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The low-energy nE~ interactions determine, in part, the role of the strange quark in dense matter, 
such as that found in astrophysical environments. The scattering phase shifts for this system 
are obtained from a numerical evaluation of the QCD path integral using the technique of lattice 
QCD. Our calculations, performed at a pion mass of ~ 389 MeV in two large lattice volumes, 
and at one lattice spacing, are extrapolated to the physical pion mass using effective field theory. 
The interactions determined from QCD are consistent with those extracted from hyperon-nucleon 
experimental data within uncertainties, and strengthen theoretical arguments that the strange quark 
is a crucial component of dense nuclear matter. 



The interactions between hyperons and nucleons are im- 
portant for understanding the composition of dense nu- 
clear matter. In high-density baryonic systems, the large 
values of the Fermi energies may make it energetically 
advantageous for some of the nucleons to transform into 
hyperons via the weak interactions, with the increase in 
rest mass being more than compensated for by the de- 
crease in combined Fermi energy of the baryon-lepton 
system. This is speculated to occur in the interior of neu- 
tron stars, but a quantitative understanding of this phe- 
nomenon depends on knowledge of the hyperon-nucleon 
(YN) interactions in the medium. In this letter we use 
nS~ scattering phase shifts in the So and 3 S*i spin- 
channels calculated with Lattice QCD (LQCD) to quan- 
tify the energy shift of the £~ hyperon in dense neu- 
tron matter, as might occur in the interior of a neutron 
star. Our results strongly suggest an important role for 
strangeness in such environments. 

Precise nucleon-nucleon (NN) interactions constrained 
by experiment and chiral symmetry, together with nu- 
merically small but important three-nucleon interactions, 
have served as input to refined many-body techniques for 
studying the structure of nuclei, such as Green-function 
Monte-Carlo [J, the No-Core Shell Model 0, or lattice 
effective field theory [3], which have led to remarkably 
successful calculations of the ground states and excited 
states of light nuclei, with atomic number A < 14. By 
contrast, the YN potentials, which are essential for a first- 



principles understanding of the hypernuclei and dense 
matter, are only very-approximately known. Therefore, 
gaining a quantitative understanding of YN interactions 
- on a par with knowledge of the NN interactions - 
through experimental and LQCD methods, is a funda- 
mental goal of nuclear science. 

Existing experimental information about the YN in- 
teraction comes from the study of hypernuclei [H [5] , the 
analysis of associated A-kaon and E-kaon production in 
NN collisions near threshold [6T-HT]. hadronic atoms |12) . 
and from charge-exchange production of hyperons in 
emulsions and pixelated scintillation devices [T3]. There 
are only a small set of cross-section measurements of the 
YN processes, and not surprisingly, the extracted scatter- 
ing parameters are not accurately known. The potentials 
developed by the Nijmegen [H [H] and Jiilich P^HTS] 
groups are just two examples of phenomenological mod- 
els based on meson exchange, but the couplings in such 
models are not completely determined by the NN interac- 
tion and are instead obtained by a fit to the available YN 
data. In Refs. [131 US] , for example, six different models 
are constructed, each describing the available YN cross- 
section data equally well, but predicting different values 
for the phase shifts. Effective field theory (EFT) de- 
scriptions have also been developed [T9l - f23] and have the 
advantage of being model independent. 

In the absence of precise experimental measurements, 
LQCD calculations can be used to constrain the YN in- 
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teractions. Several years ago, the NPLQCD Collabora- 
tion performed the first n/ = 2 + 1 LQCD calculations 
of YN interactions [24] (and NN interactions [25]) at un- 
physical pion masses. Quenched and dynamical calcula- 
tions were subsequently performed by the HALQCD Col- 
laboration [25] and by NPLQCD [27]. Recent work by 
NPLQCD [2HHS3 and HALQCD [3T1[32] has shown that 
the S = — 2 H-dibaryon is bound for pion masses larger 
than those of nature, and NPLQCD [3jl] has shown that 
the same is true for the with S = —4. In this 

letter, we use the results of LQCD calculations to deter- 
mine leading-order (LO) couplings of the YN EFT (using 
Weinberg power counting |23j ) which in turn allow for a 
determination of YN interactions at the physical pion 
mass. 

In LQCD, Liischer's method [53H56] can be employed 
to extract two-particle scattering amplitudes below in- 
elastic thresholds. For a single scattering channel, the 
deviation of the energy eigenvalues of the two-hadron 
system in the lattice volume from the sum of the single- 
hadron masses is related to the scattering phase shift 
S(q). The Euclidean time behavior of LQCD correlation 
functions of the form C x (t) = (0|x(i)x^(0)|0), where % 
represents an interpolating operator with the quantum 
numbers of the one-particle or two-particle systems under 
consideration, determines the ground-state energies of 
the one-particle and two-particle systems, E 1 ' — rn^.B 
and E^ B ^ = q 2 + m 2 A + q 2 + m 2 B , respectively. The 
form of the interpolating operators and the methodology 
used for extracting the energy shift are discussed in detail 
in Ref . [37J . By computing the masses of the particles and 
the ground-state energy of the two-particle system, one 
obtains the squared momentum q 2 , which can be either 
positive or negative. For s-wave scattering below inelas- 
tic thresholds, q 2 is related to the real part of the inverse 
scattering amplitude through the eigenvalue equation [34] 
(neglecting phase shifts in I > 4 partial- waves) : 

1 ljl<A 1 
q cot (5(g) = — lim 47rA.fl) 

J W 

This relation enables an LQCD determination of the 
value of the phase shift at the momentum ^/q 2 . 

Determining the ground-state energy of a system in 
multiple lattice volumes allows for bound states to be 
distinguished from scattering states. A bound state cor- 
responds to a pole in the S matrix, and in the case of a 
single scattering channel, is signaled by cot S(q) — > +i in 
the large- volume limit. With calculations in two or more 
lattice volumes that both have q 2 < and q cot 5(q) < 
it is possible using Eq. ([!]) to perform an extrapolation 
to infinite volume to determine the binding energy of 
the bound state = 7 2 /m, where 7 is the binding 
momentum [34-36 . The range of nuclear interactions 
is determined by the pion mass, and therefore the use of 



Liischer's method requires that m^L ^> 1 to strongly sup- 
press the contributions that depend exponentially upon 
the volume, e~ m ^ L [38] . However, corrections of the form 
e~' yL , where 7 -1 is approximately the size of the bound 
state, must also be small for the infinite volume extrap- 
olation to rapidly converge. 

Our results are from calculations on two ensembles 
of nt = 2 + 1 anisotropic clover gauge-field configura- 
tions [35] HO] at a pion mass of ~ 389 MeV, a spa- 
tial lattice spacing of b s ~ 0.123(1) fm, an anisotropy 
of £ ~ 3.5, and with spatial extents of 24 and 32 lat- 
tice sites, corresponding to spatial dimensions of L ~ 3.0 
and 3.9 fm, respectively, and temporal extents of 128, and 
256 lattice sites, respectively. A detailed analysis demon- 
strates that the single-baryon masses in these lattice en- 
sembles are effectively in the infinite- volume limit [41j . 
and that exponential volume corrections can be neglected 
in this work. Liischer's method assumes that the contin- 
uum single-hadron energy-momentum relation is satisfied 
over the range of energies used in the eigenvalue equation 
in Eq. ([!]). As discussed in Refs. [25] [3D], the uncertain- 
ties in the energy-momentum relation translate to a 2% 
uncertainty in the determination of q 2 . 

We focus on 1 S and 3 S\ interactions, NT, in the 
7 = 3/2 channel, and do not consider the 7 = 1/2 NT,- 
NA coupled channels. Calculations in the 1=1/2 chan- 
nel are complicated by the proximity in energy of the 
ground and first-excited levels in the finite volume. More- 
over, while the £ is more massive than the A, the pres- 
ence of A's in dense matter does not lower the electron 
Fermi energy. In the limit of SU(3) flavor symmetry, 
the ^o-channels are in symmetric irreducible representa- 
tions of 8 ® 8 = 27 © 10 © TO © 8 © 8 © 1, and hence 
the n£~ (along with the S _ S _ , S _ S _ , nn, and £~S~) 
transforms in the 27. YN and NN scattering data along 
with the leading SU(3) breaking effects, arising from the 
light-meson and baryon masses, suggest that all of these 
channels are attractive at the physical pion mass, and 
that S~H~ and £~£~ are bound [42H44] . By contrast, 
the 3 jSi -channel of n£ _ scattering transforms in the 10 in 
the limit of SU(3) symmetry, and is therefore unrelated 
to NN interactions. Hence, this channel is quite uncer- 
tain, with disagreements among hadronic models as to 
whether the interaction is attractive or repulsive. 

The low-energy nS~ interactions can be described by 
an EFT of nucleons, hyperons and pseudoscalar mesons 
(71-, K and 77), constrained by chiral symmetry 19!.i22 ] [23 ] . 
At leading order (LO) in the expansion, the nS~ inter- 
action is given by one-meson exchange together with a 
contact operator that encodes the low-energy effect of 
short-distance interactions. As these contact operators 
are independent of the light-quark masses, at LO the 
quark-mass dependence of the interactions is dic- 

tated by the meson masses. Therefore, in each partial 
wave, a single lattice datum at a sufficiently low pion 
mass determines the coefficient of the contact opera- 
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FIG. 1: LQCD-predicted % nE" phase shift versus labo- 
ratory momentum at the physical pion mass (blue bands), 
compared with other determinations, as discussed in the text. 
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FIG. 2: LQCD-predicted 3 Si nS~ phase shift versus labo- 
ratory momentum at the physical pion mass (blue bands), 
compared with other determinations, as discussed in the text. 



tor, thereby determining the LO interaction, including 
energy-independent and local potentials, wavefunctions 
and phase shifts, at the physical pion mass. 

We find that our LQCD calculations in the 7iX~ 
channel are consistent with the SU(3) symmetry expecta- 
tions. At m n ~ 389 MeV, using a volume extrapolation 
as discussed above, we find that this channel has a bound 
state, with binding energy B = 25 ± 9.3 ± 11 MeV. The 
quality of the LQCD data in the nS" channel is com- 
parable to that of its 27-plet partner analyzed in 
detail in Ref. [3D] (see also 05]). In the EFT, the coeffi- 
cient of the LO contact operator in this channel is deter- 
mined by tuning it to reproduce the LQCD-determined 
binding energy. We find that this channel becomes un- 
bound at m^tS 300 MeV, in agreement with Ref. [46] . 
which constrained the LO contact operator using exper- 
imental data. In Fig. [TJ we show the predicted -"So nTT 
phase shift at the physical pion mass — (dark, light) blue 
bands correspond to (statistical, systematic) uncertain- 
ties — and compare with the EFT constrained by ex- 
perimental data [23], the Nijmegen NSC97f model [14], 
and the Julich '04 model [IS]. The systematic uncer- 
tainties on our predictions include those arising from the 
LQCD calculation (see 05J) as well as estimates of omit- 



ted higher-order effects in the EFT. 

The 3 Si- 3 Di nY>~ coupled channel is found to be highly 
repulsive in the s-wave at ~ 389 MeV, requiring in- 
teractions with a hard repulsive core of extended size. 
Such a core, if large enough, would violate a condition re- 
quired to use Liischer's relation, namely R <C L/2 where 
R is the range of the interaction. We have determined 
the EFT potential directly by solving the 3-dimensional 
Schrodinger equation in finite volume to reproduce the 
energy levels obtained in the LQCD calculations. The re- 
pulsive core is found to be large, and formally precludes 
the use of Liischer's relation, but both methods lead to 
phase shifts that agree within uncertainties, indicating 
that the exponential corrections to Liischer's relation are 
small. In Fig. [2] we show the predicted 3 Si nS~ phase 
shift at the physical pion mass. 

The nTi~ interactions presented here are the crucial 
ingredient in calculations that address whether £~'s ap- 
pear in dense neutron matter. As a first step, and in order 
to understand the competition between attractive and re- 
pulsive components of the interaction, we adopt a 
result due to Fumi for the energy shift due to a static 
impurity in a non-interacting Fermi system |47j : 



AE 



1 

7T/X 



dk k 



(2) 



where /i is the reduced mass in the riT,~~ system. Us- 
ing our LQCD determinations of the phase shifts, and 
allowing for a 30% theoretical uncertainty, the resulting 
energy shift and uncertainty band is shown in Fig. [3] At 
neutron number density p n ~ 0.4 fnV 3 , which may be 
found in the interior of neutron stars, the neutron chem- 
ical potential is \i n ~ Mn + 150 MeV due to neutron- 
neutron interactions, and the electron chemical potential, 
fjt e - ~ 200 MeV 08]. Therefore fx n + fx e - ~ 1290 MeV, 
and consequently, if /x E - = M E + AE < 1290 MeV, that 
is, AE < 100 MeV, then the E~, and hence the strange 
quark, will play a role in the dense medium. We find 
using Fumi's theorem that AE = 46 ± 13 ± 24 MeV at 
p n = 0.4 fnV 3 . Corrections due to correlations among 
neutrons are difficult to estimate and will require many- 
body calculations which are beyond the scope of this 
study. Despite this caveat, the results shown in Fig. [3] 
indicate that the repulsion in the n£ - system is inad- 
equate to exclude the presence of S~'s in neutron star 
matter, a conclusion that is consistent with previous phc- 
nomenological modeling (for a review, see Ref. 09]). 

In this letter, we have presented the first LQCD predic- 
tions for hypernuclear physics, the ISo and 3 Si n£ _ scat- 
tering phase shifts shown in Fig. [T] and Fig. [2] While the 
LQCD calculations have been performed at a single lat- 
tice spacing, lattice-spacing artifacts are expected to be 
smaller than the other systematic uncertainties. We an- 
ticipate systematically refining the analysis presented in 
this letter as greater computing resources become avail- 
able. The n£ _ interaction is critical in determining the 



4 




FIG. 3: The energy shift versus neutron density of a single 
£~ in a non-interacting Fermi gas of neutrons as determined 
from Fumi's theorem in Eq. (J2|. The inner (outer) band en- 
compasses statistical (systematic) uncertainties. 



relevance of hyperons in dense neutron matter, and we 
have used the LQCD predictions of the phase shifts to es- 
timate the £~ energy shift in the medium. Our calcula- 
tion suggests that hyperons are important degrees of free- 
dom in dense matter, consistent with expectations based 
upon the available experimental data and hadronic mod- 
eling. It is important that more sophisticated many-body 
techniques be combined with the interactions determined 
in this work to obtain a more precise determination of the 
energy shift of the E - in medium. This will refine the 
prediction for the role of strange quarks in astrophysi- 
cal environments, and, in particular, will quantitatively 
address questions posed by the recent observation of a 
1.9 Mq neutron star [50] . 
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